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General Process for Developing and Implementing 
a Regulated Biochemical Reaction Network Model 



START ^) 



110 

Construct a biochemical reaction network 
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Mathematically describe the reaction network 



Determine the regulatory requirements for reactions in the 
network to take place 
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140 

Mathematically describe the regulatory controls 



150 

Formulate problem by specifying parameters and initial 
reaction constraints to the model 



160 

Determine if regulation limits any reactions 
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Refine reaction constraints 
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Simulate the behavior of the reaction network subject to all 
the constraints 
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A) 



Example Biochemical Reaction Network 




B) 



Example Regulatory Structure and Requirements 



C) Simulated Reaction Network (without regulation) 




D) 



Simulated Reaction Network (with regulation) 
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Metabolic Model 
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Regulatory Model 
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•Activity constraints set for rxn Itoith 

(lower bound = 0, upper bound = INF or #) 

•Inactivity constraints for rxn stolch 
(lower bound = 0, upper bound = 0) 



Integration of Stoichiometric model and 
Logical model achieved through regulatory 
restraints (logic values of reaction 
processes) which are used to refine 
appropriate reaction constraints in the 
model: 

If rxn Logic = 1 then use Activity constraints 
If rxn^ . = 0 then use Inactivity constraints 



Logic functions 

a, = (activator/inhibtor) • TF 

a 2 =l 

c, = TF* • pr, • genel * gene2 
C 2 = P>3 * gene3 

•l = M gencl 

12 = M 8 cnc2 

1 3 = M gcne3 

pi = Pgencl * Pgene2 " ^gene3 

rxn Lx) g ic = Protein • Cofactor • Substrate, • Substrate 2 

Time delays can be specified for the switching 
of each memorization variable after a 
triggering change in the associated function 
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200 

A Time-Dependent Implementation of a Regulated 
Biochemical Reaction Network Model 
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210 

Determine the time steps to be used in the transient 
calculations for the simulation 



220 

Formulate problem by specifying the intial time conditions 
and parameters 



230 

Determine status of regulatory variables for the time point 
under consideration 



Formulate problem by specifying the conditions and 
parameters for the next time point 



240 

Refine the constraints on the biochemical reactions based 
on the status of their regulation 



Use the flux distribution to determine the new set of intial 
conditions for the next time point 
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Solve for a steady-state flux distribution in the reaction 
network subject to the constraints 
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300 

Process for Developing Genome Scale Regulated 
Models of Metabolism 
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Determine the location of all the open reading frames in the 
genome 
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Annotate the open reading frames with functional 
assignments 
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Determine genes involved in cellular metabolism and/or 
metabolic regulation 



351 

Gather biochemical information on all reactions performed 
by the genes£ene products of the metabolic genoype 
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Group metabolic and regulatory genes into transcriptional 
units based on sequence analysis 
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352 

Gather additional biochemical information on any reactions 
known to occur biochemically but not present in metabolic 
genotype 
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Locate regulatory binding sites on transcriptional units by 
sequence or tnctional homology 










353 

Generate an organism specific reaction index and 
stoichiometric matrix 




363 

Gather biochemical information on all regulation of 
transcriptional units 
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354 

Specific default reaction constraint conditions based on 
thermodynamics and reaction kinetics where available 




364 

Define the organism specific regulatory structure using 
Boolean logic representation 
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Determine organism specific biochemical and physiological 

data on biomass composition, uptake rates and 
maintenance requirements 



380 

Formulate the general linear programming problem 
combined with regulatory constraints representing the 
organism specific model 
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